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INTRODUCTION 

Fatigue  fractures  occur  at  points  of  stress  concentration, 
which  are  normally  regions  of  complex  stress  distribution  with 
large  gradients.   The  stress  field  is  geometry  dependent  and 
will  only  be  reproduced  in  a  laboratory  sample  if  geometric 
similitude  is  preserved.   Matching  sample  to  part  must  also 
include  similitude  in  the  two  major  phases  of  fatigue;  initia- 
tion and  propagation  and  not  merely  a  matching  of  stress  con- 
centration factors.   In  fact  if  the  fatigue  process  were  well 
understood,  it  would  seem  that  a  knowledge  of  the  local 
conditions,  coupled  with  material  property  information,  would 
be  all  that  is  required  to  predict  fatigue  life. 

This  has  not  been  possible  so  far,  but  this  study  has  been 
undertaken  to  see  how  much  success  can  be  achieved  toward  this 
end  via  a  detailed  study  of  conditions  surrounding  the  notch 
tip.   To  have  any  chance  of  success  whatsoever,  the  local  stress 
and  strain  conditions  at  the  notch  must  be  determined  accu- 
rately.  One  of  the  means  available  for  accomplishing  this 
is  experiment,  and  this  approach  was  used  successfully. 

SIMILITUDE 
An  evaluation  of  how  well  one  notch  specimen  models 
another  can  be  studied  via  similitude  theory.   The  strains  in 
the  vicinity  of  a  notch  depend  upon  the  following  variables: 

ti  -  \(  p  w,  t,yO,  A,  e,o)  (1) 


where  £•  =  Strain  components 

P  =  Applied  load 

W  =  Sheet  width 

t  =  Sheet  thickness 

P  =  Radius  of  curvature  of  notch 

a  =  Depth  of  notch 

£.  =  Tensile  Modulus 

V  =  Poisson's  Ratio 

* 

By  dimensional  analysis,  the  functional  relationship  can  be 
written  in  terms  of  groupings  of  variables  that  are  dimen- 
sionless . 


where 

P 


(T  = 


Wt 


From  similitude  theory,  if  in  two  different  geometries 


0,  =  ^v 


Then 


GEOMETRIC  SIMILARITY 

If  the  material  of  the  part  and  the  model  are  the  same, 
the  same  strains  are  produced  only  if  they  are  geometrically 
similar.   It  is  not  sufficient  to  match  only  the  stress  at 
the  notch  tip,  but  complete  similarity  must  be  maintained. 
If  the  complete  strain  fields  are  not  the  same,  the  plastic 
zones  will  be  different,  and  the  neighborhood  where  fatigue 
action  is  taking  place  has  different  conditions.   Thus, 
similarity  must  be  preserved  when  employing  experimental 
fatigue  results  from  one  specimen  to  predict  fatigue  failure 
in  another.   This  is  not  generally  done  in  current  practice. 
Since  it  is  not  convenient,  nor  even  practical,  to  run  tests 
on  a  geometrically  similar  specimen  each  time  a  part  is  an- 
alyzed for  fatigue,  a  different  approach  is  attempted  to 
develop  a  more  general  theory  where  knowledge  of  local 
conditions  at  the  notch  tip  is  sufficient,  along  with 
material  properties,  to  predict  fatigue  in  any  specimen. 

SIMILITUDE  IN  UNIAXIAL  TENSION 

Because  much  of  the  fatigue  data  available  has  been 
collected  on  uniaxial  tension  specimens  and  because  some- 
thing resembling  uniaxial  tension  prevails  at  notch  tips, 
it  is  reasonable  to  pursue  this  similarity,  which  is  shown 
in  Figure  1. 
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Figure  1.   Uniaxial  Model  of  Stress  Condition  at 
a  Notch  Tip. 


For  an  elastic  material 


tL--  k^,^,^s°) 


o: 


If  the  same  material  is  used  for  the  model  and  the  specimen, 
the  strain  field  will  be  modeled  to  the  degree  that  the 
stress  field  is  modeled.   For  instance,  in  the  region  of 


the  notch  tip  in  plane  stress, 


U  --  \  [o;  -  ^<ry]  (4) 


* 


To  keep  the  departure  from  a  uniaxial  stress  state  less  than 
3%,  which  would  be  in  the  range  of  experimental  measurement 
error,  the  transverse  stress  must  be  small. 


(j;    <:     o.o*>dG) 


(5; 


For  a  typical  value  of  Poisson's  Ratio  for  aluminum,  V  =  0.3  3 
and 

d      <         O.0\    Cy 

In    the   y   direction, 


6) 


t,.if<r,.o<r.] 


(7) 


To  maintain  the  same  3%  limit, 

0<u      i   0.03  G} 


or 


Also  in  the  Z  direction, 

(7*  <    ^-^  °) 


(8) 


(9) 


The  most  demanding  of  these  is  the  first.   In  order  to 
model  the  behavior  at  the  notch  tip  as  a  uniaxial  specimen, 
only  regions  may  be  considered  where  the  radial  stress  is  one 
per  cent,  or  less,  of  the  hoop  stress. 

For  a  circular  hole  in  an  infinite  sheet,  the  location 
where   GV  =  O.OI  G"©   is  at  f*  -  I.OIOI  CL  .    For  a  two 
inch  diameter  hole,  the  distance  from  the  edge  of  the  hole 
where  the  error  exceeds  the  stated  allowable  is  0.01  in. 
This  precludes  the  use  of  a  strain  gage  to  make  local  measure- 
ments.  If  specimen  sizes  were  enlarged  ten  times  until  the 
uniaxial  tensile  region  were  large  enough  to  accommodate  a 
1/8"  strain  gage,  the  hole  would  be  twenty  five  inches  in 
diameter  and  the  specimen  would  be  four  to  six  feet  wide  and 
ten  to  twelve  feet  long.   To  circumvent  this  difficulty, 
photoelasticity  was  used  as  the  experimental  method  for 
measuring  local  strain  conditions  at  the  notch.   There  are 
similar  size  problems,  but  they  are  not  quite  as  severe  as 
for  strain  gages,  plus  the  bonus  of  obtaining  the  entire 
strain  field  at  the  notch  tip  made  photoelasticity  attractive 

PHOTCELASTIC  MATERIAL  PROPERTIES 

SELECTION 

70  75-T6    aluminum  was    selected   for    the   study   of   notch 
fatigue   since   it  is    a  popular   alloy  of    the    aircraft  industry 
and  one   for  which    there   is    a   sizeable    fatigue    data  base. 


To   coat    7075-T6   with    a  birefringent  material    that  may  be 
loaded   into    the   plastic    zone,    Photolastic* s   PS-1  material  was 
selected.      It  has    a  maximum   strain   capability   of    10%    and   is    a 
good   all    around  photoelastic  material. 

Mil  Handbook.    5    gives    typical   data    for    7075-T6    sheets 
of    0.080    in    thickness.      To   extend   into   the   yield   range,    strains 
of    10,000   ys   must  be   measured  with   excursions    up    to    20,000    ys 
probably   anticipated.      From  equation    (A-l)    of   Appendix  A,    it 
is    readily   seen    that 

Ft    —    _  do) 

If  something  like  five  fringes  are  desired  at  the  10,000  ys 
level,  the  F«  value  needed  is  0.002.   This  would  allow  10 
fringes  at  the  maximum  strain  level  expected. 

From  equation  (A-4) ,  the  material  fringe  value  is  given 
bv 

For  PS-1,  the  company  supplied  a  value  of  K  =  0.15,  and  for 

white  light  the  wave  length  is  r  =  22,7  x  10~   in;  thus 

-4 
the  material  fringe  value  is  1.5  x  10   .   Using  equation  (A-l) , 

the  thickness  of  the  coating  desired  is 

4 


The  coating  material  PS-lC  was  purchased  in  sheets  nominally 

.040  in  thick.   The  bonding  agent  used  was  that  suggested  by 

Photolastic.   At  first  PC-1  was  used,  but  for  higher 

elongations  it  was  necessary  to  go  to  PC-8,  even  though  the 

first  cement  was  supposed  to  have  the  elongation  capability 

required. 

CHARACTERIZATION    OF    PS-lC    PH OTOE LAS TIC    MATERIAL 

A   careful   characterization   of   PS-lC  was    conducted  by 
Stenstrom    (1)    using    two   uniaxial    tensile    test   specimens 
prepared   from   a   single   sheet  of    the   photoelastic  material. 
The    first  specimen  was    1.0    inch  wide    and  was    used   for   thick- 
ness  measurement  and   longitudinal   strain    data.      The    second 
2.0    inch  wide    specimen  was    used   for    transverse    strain    data 
and   for   measuring    the   photoelastic   coating    stress    fringe 
value,    F_  . 
Thickness    Determination 

One    longitudinal  edge   of    the    1.0    inch  wide    specimen  was 
lightly   sanded   to    create   a   smooth,    uniform  cross-section   of 
the   photoelastic  material.      It  was    then  mounted    such    that   a 
traveling  microscope    could  be    used   to  measure    thickness.      With 
the    contrast  in    appearance   between   the    reflective   coating    and 
the   base   material,    it  was   possible    to   accurately   determine 
the    thickness    of   each.      For    the    specimen   used,    the   base 
thickness  was   measured   to  be    0.0382    inches    and   the    reflective 


coating  to  be  0.0051  inches.   These  numbers  were  obtained  by 
measuring  eight  separate  points  along  the  edge  and  averaging 
the  results.   The  total  average  thickness  of  0.0  433  inches 
fell  within  the  range  of  the  manufacturer's  measurement  of 
0.044±0.002  inches.   It  was  therefore  concluded  that  the 
labeled  thickness  included  both  the  base  and  the  reflective 
coating  after  two  similar  specimens  of  different  marked 
thicknesses  were  checked.   For  calculation  purposes,  the 
thickness  of  the  coating  was  taken  to  be  0.0  382  inches  since 
the  path  of  polarized  light  does  not  penetrate  the  reflective 
coating  of  the  phctoelastic  material.   Inherent  in  using 
this  thickness  is  the  assumption  that  the  reflective  coating 
does  not  contribute  to  the  reinforcement  of  the  base  material 
Young ' s  Modulus 

Two  lines  were  lightly  etched  on  the  surface  of  the  1.0 
inch  wide  specimen  in  the  horizontal  direction  in  order  to 
accommodate  a  Tuckerman  Optical  Strain  Gage.   The  specimen 
was  then  loaded  in  tension.   The  readings  were  taken  with  a 
autocoliimater  and  recorded.   Using  the  American  Instrument 
Company  calibration  factors,  the  strain  was  calculated  with 
the  following  formula  taken  from  Ref.  2. 


b  =    - — d2) 


loooe 


where  S  =  Strain  in  microinches  per  inch 

F  =  Lozenge  calibration  factor,  =  1.002 

A  =  Auto-collimater  calibration  factor,  =  1.00  4 

L  =  Lozenge  length,  2.0  inches 

R  =  Net  reading 

e   =   Gage    length,    1.0    inch 
With    the   measured   strain,    the  measured    load   at  each   point 
and  the  known   cross-sectional   area,    Young's   Modulus,    E,    was 
determined   from   the   slope    of    the    stress-strain   curve    shown   in 
Figure    1,    using    a   linear    regression   analysis. 

<Jz     3S7,0Ti.»2  £,  -  4.*w>i  (i3) 

fc_;      0.9966 

where   R  =   correlation   coefficient. 

Dropping   the   non-zero   intercept,    Young's   modulus   was    taken 

to  be 

£cr     5S1,079      psi.  (14) 

Photolastic  Inc.  information  lists  this  modulus  to  be  nomin- 
ally 360,000  psi.   Using  thin  specimens  for  this  measure- 
ment requires  very  accurate  thickness  measurements  because 
E  is  inversely  proportional  to  thickness.   Small  errors  in 
thickness  can  turn  out  to  produce  large  percentage  errors  in 
area  and  modulus. 
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PS-1C    CHARACTERIZATION 
STRESS-STRAIN    CURVE 


1500 


1000    - 


Stress, 

a   psi 


500    _ 


Strain,    e,    10"      in/in 


Figure    1 
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Poisson's    Ratio 

Using    the    2.0    inch  wide   PS-lC   specimen,    two    lines   v/ere 
lightly   etched   into   the    surface    in    the    longitudinal   direction 
one    inch   apart.      With    the   Tuckerman   Optical  Strain   Gage 
mounted   in   the   transverse    direction   as    shown    in   Figure    2, 
the    specimen  was    loaded  in    tension.      The    data  points    from 
the    test  were    used   in    a    linear   regression   analysis   with    the 
following    results: 

0-.  -941^87.02.  £t  ~  l.ooioe  (15) 

R,=    0.99991 

Poisson's    ratio  was    then   determined  by   dropping    the   non-zero 
intercept   and   dividing    the   negative   of    the    slope   by  E    as 
indicated  in   Figure    1.      The   measured  Poisson's    ratio   is 

0    =    0.  1195 

which  compares  with  a  nominal  value  of  0.3  8  supplied  by 
Photolastic  Inc.   These  tests  were  performed  in  accordance 
with  ASTM  standards  for  uniaxial  properties  (3) . 
Coating  Fringe  Value  in  Terms  of  Stress,  Fa 

Prior  to  etching  lines  on  the  2.0  inch  wide  specimen,  it 
was  loaded  in  tension  in  the  elastic  region,  and  the  fringe 
number,  N,  was  determined  using  a  reflection  polariscope. 
With  several  load  points  and  known  dimensions,  the  stress 


12 
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was  calculated  and  plotted  against  the  fringe  number.  Using 
the  data  from  Table  I,  a  linear  regression  was  made  with  the 
data  which  yielded 

(T-    4-Slo.L!*  N  ■+  l.9£14  (16) 

Dropping  the  non-zero  intercept,  the  stress  fringe  value  was 
taken   to  be 


F^r     4-SC.tl        P^l/fv^gCL  (17) 

Coating   Fringe   Value   in   Terms    of   Strain,    F£ 

With    the    tensile  modulus    and  the   Poisson's    Ratio   of 
the   coating   material   already   determined,    the    strain   coating 
fringe   value   can  be   calculated  readily.      From  Appendix  A, 


CHARACTERIZATION   OF    70  75-T6 
The   elastic   and  plastic  material  properties   of    the 
aluminum  panels  were    established  by  KAISER    (4)    using    tensile 
tests   of    uniaxial   specimens   made    from   the    same   mill    run   as 
the    sheets    to   be    tested.      The   specimens   were   manufactured 
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TABLE     I 


PS-1C      F         DETERMINATION 


0 

psi 

Compensator 
Reading 

N 

0.0 

0.0 

0.0 

202.37 

23.0 

0.489 

408.42 

39.0 

0.330 

610.26 

63.0 

1.340 

785.26 

79.5 

1.691 

990.53 

102.0 

2.170 

189.21 

123.0 

2.617 

LINEAR    REGRESSION: 

o    =    456.6661N    +    1.92740 
R   =    .99932 
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and   tested   according   to  ASTM  standards  E8    and  B5  57.     (3)       The 
specimens  were    cut   transverse   to    the   rolling   direction,    and 
the   tests  were  performed   in   a    10  0  Kip  MTS   machine.      The 
duration   of    the   tests  was    three  minutes.      MICRO-MEASUREMENTS, 
EA-13-125AD-120 ,    precision    strain  gages   with   a   temperature 
compensated  bridge   circuit  were   used   on   all   specimens.      Trans- 
verse gage   sensitivity  errors   were   corrected   according    to   the 
manufacturer's    recommendations    (5).      Critical   cross-section 
measurements  were  made  with   a  micrometer. 
YOUNG'S    MODULUS     (E) 

Tests   were   conducted  using    the    specimen    shown   in 
Fig.    3   on   the   MTS    test  machine   with    10,000    lb   INSTRON   grips, 
which   gripped   the   specimen   evenly.      The    results    of    testing 
three    specimens    are   shown   in   Tables   II    to   IV  where    all   strain 
values   have   been   corrected   for   transverse    sensitivity. 
Linear   regression   of    the   first  eight  points    in    the   elastic 
range   of    all   three   specimens   gave 

G"  =     i%L  %U  +•   Lon7  *  lo1  £•  (18) 


2 

with   a   correlation   coefficient   of    R     =   0.9997.      Ignoring 


the   non-zero   intercept   as    specified  by   ASTM 


£  r       I  0.  I  Z.  X  lOfe    p*L  (19) 
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r~\    / 


4.0' 


8.25" 


t_; 


0.5 


0.75 


AXIAL     STRAIN 
GAGE 


0.75     R 


ThicHness  =  0.0795 


FIGURE  3     Characterization     Specimen 
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TABLE    II.       MTS    SPECIMEN    A   TEST    RESULTS 


Cross-section 
Load,    lbs. 


0.03975    in2 


-6 


Strain,    £,10        in/in 


256 
503 
750 
1005 
1255 
1505 
1778 
2003 
2252 
2508 
2755 
2905 
2984 
3037 


615 
1204 
1801 
2423 
3042 
3665 
4352 
4925 
5587 
6355 
7365 
8230 
9045 
10150 
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TABLE    III.       MTS    SPECIMEN    B    TEST    RESULTS 

2 

Cross-section   =   0.03975    in 

Load,    lbs. Strain,    z,    10        in/in 

500  1250 

750  1860 

1000  2450 

1250  3080 

1500  3700 

1750  4320 

2000  4950 

2250  5630 

2500  6380 

2750  7340 

2900  8200 

3000  9250 

3050  10200 

3100  11550 

3125  12800 
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TABLE    IV.       MTS    SPECIMEN    C   TEST    RESULTS 

2 
Cross-section   =   0.0  39  75    in 

Load,    lbs. Strain,    s,    10        in/in 

272  663 

503  1220 

76-2  1850 

1008  2451 

1231  3005 

1506  J690 

1755  4318 

2000  4940 

2255  5607 

2503  6325 

2750  7230 

2900  8120 

3000  9200 

3060  10500 

3100  11500 

3125  12500 
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POISSON'S    RATIO    (v) 

Tests   were   conducted   using    the    specimen   shown    in 
Fig.    4   on    the    RIEHLE    test  machine   with    10,000    lb    RIEHLE 
grips.      The    results   were    tabulated   in   Table   V.      The    trans- 
verse   strain   measurements   were    corrected  by    the   manufacturer's 
data  on    transverse    sensitivity.      Both    strains   were    read   at 
the    same    load   and  plotted  versus    one    another.      A    linear    re- 
gression was   made   en    the    data   resulting   in: 


f       -       O  OOOO    -    O.hZSt*   £-, 


(20) 


with  a  correlation  coefficient  of  R  =  0.99998.   Disregarding 
the  constant  factor  as  insignificant  to  four  places 


$  z      0.  ^ZSl* 


(21) 


YIELD   STRESS    AND   STRAIN    HARDENING    MODULUS 

The   values   of   yield   stress    and   strain  hardening   modulus 
were    determined  graphically   using    the    data    from   Tables   III 
and  IV.      Plastic    region   data   in   Table   II    is    not   reliable 
because   of   excessive    creep    encountered   during    that   test. 
0.2%    offset  yield   stress,    ov   =    76,000    psi 
strain  hardening   modulus,    E^.    =    566,000    psi 
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FIGURE     4 

2    GAGE    SPECIMEN 


11.72 


1.5" 


AXIAL    STRAIN 
GAGE 


TRANSVERSE 
STRAIN      GAGE 


thickness  =    0.0800 
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TABLE    V.       REIHLE    SPECIMEN    TEST    RESULTS 

2 

Cross-section   =   0.12    in 

LOAD 
lbs. 

500 
1000 
1500 
2000 
2500 
3000 
3500 
4000 
4500 
5000 
5500 
6000 
6500 
7000 
7500 
8000 
8500 
9000 
9500 


STRAIN, ex 

10"6in/in 

345 

720 

1114 

1501 

1895 

2291 

2705 

3095 

3515 

3912 

4332 

4750 

5188 

5595 

6075 

6649 

7385 

8663 

12245 

STRAIN, £2 

10"6in/in 

-115 

-240 

-36  4 

-492 

-614 

-745 

-871 

-993 

-1126 

-1255 

-1384 

-1520 

-16  4  7 

-1784 

-1928 

-2108 

-2375 

-2888 

-4435 

23 


The   graphical   fit  of    these   values    to    the    test    data   can  be 
seen   in   Fig.    5 . 
RAMBERG-OSGOOD    COEFFICIENTS 

The    Ramberg-Osgood   equation   for  elastic-plastic 
stress-strain    characterization   is    given  by: 


t-Z*Mi 


(22) 


where: 

£  =  strain 

(T  =  stress 

E.  =  Young's   modulus. 
The   8    and   n  coefficients   were    determined  graphically   from 
the    data  of    Table   IV,    by    the   method  given   by   Ramberg    and 
Osgood    (6)  . 


A 


1     L        <Tl6    J 


(23) 


r\  --    I 


^  (7.) 


3  n  cr.«£  J 


(24) 


Lines  corresponding  to  slopes  of  0 . 7E  and  0 . 85E  were 
laid  out  on  the  X-Y  recorder  plots  of  load  versus  longitudinal 
strain.   Load  values  where  the  rays  intersected  the  load- 


24 


92 


80 


FIGURE  5 
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£=  10.12  xlO  "  psi 
j3-   1.479  x  I043 
n*  21.58 

a  -  SPECIMEN  a 

a  -  SPECIMEN  C 

SOLID  LINE  -  RRMBERG-QSSCOD 
CURVE 


%  offset 


a 


20 


5TRRIN  ,  6  ,    10°  in/in 
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strain   curve   were 


p     =    3060   lbs  PISE  =   2  950    lb* 

«1£ 


2 
The  cross-sectional  area  from  Fig.  3  is  0.0  39  75  in  .   This 


gave 


csv  77^4a4  psl       ^™~-  14'114  P" 


The  Ramberg-Osgood  Parameters  are  given  by 

A  =    1.479   x-  I04i 

The    stress-strain    law   for   uniaxial   tension    throughout    the 
entire   range    is 

-8^-  ~    .^       ,  -103    _    Z1.S8 

£;9.88lxio    (T    +.9949X10      <j-  (25) 

This  is  shown  plotted  in  Fig.  5  along  with  data  points  from 
two  of  the  specimens. 

REINFORCEMENT 
When  specimens  are  made  from  thin  sheets  of  material 
to  duplicate  wing  skins,  the  coating  material  can  provide 
some  reinforcement  to  the  specimen  and  influence  the  stress 
values.   This  is  usually  small,  but  since  it  is  calculable 
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by  classical  developments,  allowance  can  be  readily  made 
for  it. 

Consider  an  infinitesimal  element  in  a  biaxial  stress 
field  oriented  along  principal  stress  directions  as  illus- 
trated in  Fig.  6.   If  the  total  loads  applied  to  the 
element  are  P,  and  P.,,  they  can  be  equated  to  the  loads  in 
the  coating  and  the  specimen, 

P,  =  P|»  +  P^  P*  ""  P"  *  P"       (2«) 

If  there  were  no  coating,  the  stress  produced  by  P.  and  P~ 
would  be 


P,  =  <Hu^^k  pl*-  °lH.Wsa.K, 


(27) 


The  stresses  in  the  component  parts  of  the  coated  specimen 
are  given  by 

P15>  =<r,4ws<ix->.  P!b '-  ^  Hs i* <       (28) 

Pit  -.  cr,c wt  a*2  plt  --  <*"«■ K' *x i     <29> 

Equating  the  loads  in  the  coated  and  the  uncoated  elements, 

G;Kl»,4xi*   C1SM*«-  +  G~lcKoaKt  (30) 
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II 


p       p 

r»C        rIJ 


FIGURE  6.     Reinforcement  Model 
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Simplifying    terms, 


cr,«  --  o"»  *   ^  ^ 


2.VL 


(32) 


(33) 


Subtracting  the  two  equations 


C^-(T^-C^-^)^(^,-^) 


(34) 


From  the  plane  stress  constitutive  equations, 


<TlS. 


i-o 


r(M0tl)       cr,.^v(^t, 


(35) 


equation    (34)    can  be    transformed   into  expressions    for   strain 


UrO^ 


Mx~    Cj/w. 


l  +  Oj  v  Ks   \+^c 


(36) 


Recalling  the  fundamental  characteristic  of  the  photoelastic 
coating  method. 

L^   5  £.c  £*s  -  ^tc  (37) 


t$ 


The  reinforcement  equation  simplifies  to 

I  "VL.    ^^  T~~   T — '  7 C""\ 


(38) 
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The  quantities  in  the  bracket  have  all  been  measured;  there- 
fore all  photoelastic  measurements  must  be  multiplied  by 
this  factor  to  obtain  uncoated  specimen  response 


i^-t^  '-   l.©ici  Pt^ 


39) 


POISSON'S  RATIO  MISMATCH 
As  mentioned  previously  the  notch  tip  may  be  considered, 
under  certain  restrictions,  as  a  uniaxial  specimen.   As 
the  tip  is  strained,  the  lateral  contraction  is  different  in 
the  specimen  and  in  the  coating.   At  the  interface  of  the 
coating  and  specimen, 


Ut^   *«  =  -^.*is=^s*i 


(40) 


at  the  free  surface  of  the  coating 

ttc  =  "^c  £10  (41) 

Thus  there  is  a  transverse  strain  gradient  through  the  thick- 
ness of  the  coating.   The  birefringence  will  exhibit  the 
average  strain  readings,  which  will  lie  somewhere  between  the 
values  at  the  interface  and  the  values  at  the  surface.   At 
the  interface, 
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t^'iu.   =    ^c(lA)-   FtN 


(42 


At  the  free  surface, 


fcu."  * 


LC- 


-   tic  0  +  ^0  ~  ^ 


N 


(43) 


The  actual  fringe  value,  N,  will  be  an  average  bounded  by 
these  two, 


(44) 


Experiments  with  uniaxial  tensile  specimens  have  shown  (7) 
that  the  strain  variation  through  the  thickness  are  only 
present  near  free  boundaries  in  the  plane  of  the  specimen, 
such  as  the  notch  surface.   This  arises  because  of  the  edge 


effect. 


f 


T^TX 


£iC=  "^C^iC 


tu'-dttxt. 


y 


CtlOSS-SSCTiOrJ      Ofi       UWi**lAl-    SPSCiMeM 


w 


Describing  the  displacement  field  sketched  above,  we  assume 
it  to  be  of  the  following  form: 
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U--U.(x;i)         V=V*(y)        lOzVJ^l)        (45) 


w 

At  x  =  /2  there  exists  a  contour  for  u  which  describes  a 


curve  in  z 


u(£,i)  =  gW  <46) 


2      2 
For  instance  g(z)  =  (z   -  h   )  yields  a  good  solution; 

however  for  now  we  take  a  more  general  approach  and  only 

assume  that  g(-h  )  =0. 

c 

This  contour  is  continuously  modified  as  one  proceeds  to  the 
inner  regions.   This  will  be  represented  by 


Since  the  displacement  in  the  coating  is  not  zero  at  the 
specimen  boundary  but  rather  produces  a  constant  strain  equal 
to  that  of  the  specimen,  another  term  must  be  added. 

uU,t)  =  +WgU)  +  c'x  <48» 

The  other  displacements  are  chosen  to  give  constant  strains 

v(y)  -  <uy  (49) 

-10  Cl)    :  ^^  (50) 


32 


The    strains   become 


6yi 


-   C: 


51     ^7  " 


(51) 


The  stresses  become 


(u-. 


+  oc 


f'g  +c.]  <->c[4'g  *  c,  *ct+c5^ 


C,  - 


Ec 


»    l+0c 


<r»  *  - 


£c 


4      ItOc 


(52 
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?«ls   G  +  S' 


n  --    o 


7-  -    o 


The  general  form  of  the  stress  expressions  is: 


OT.     ~    FWfi(0    *    Cor\t>ta*t  (53) 


Using  this  general  form  as  a  guide,  a  more  general  solution 
is  sought  using  the  field  equations.   Assume   xXy  =  txz  =  0 

ft 

and  t-  =0.   The  Equilibrium  Equations  simplify  to 

5*  3e 


These  equations  can  be  solved  with  a  stress  function  by 
prescribing  that 

(T  -    £l  ^*  =  -23L     ^i-  Zl    (55) 
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The    Compatibility  Equations    reduce    to 


*  l+Vc     2)*> 


Vo  (56 

Adding    the    first   three   compatibility   equations   gives 


lince  7l(Ty  -O 


VC6     =      O  (57) 


Then 


^7l(cTx  +  crs)  =  o  ,58) 

Substituting    the    stress    functions    into    (5  8) ,    we    obtain    the 
governing  equation   for   it. 

V4£     =   ©  (59) 

Using  the  biharmonic  theorem, 
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(60) 


where 


(61) 


5  ^tfi  +-  / 

V2^  =  O  and         VZ/  =0 

Using   separation   of  variables    for    $£      and     ^ 

the    stress    iunction    can  be    represented   as 

£(x,*)^   if(x)gte)  +■  KMk(2) 

The    stress    of  principal   concern  here    is    <T^  . 

a;  ^  -fW  dl[2g(^)J   +  kw  k*(i)  (63) 

To  make  it  match  the  general  form  of  equation  (5  3) ,  let 


(62) 


a  i1 


(64) 


g(z)    is    found  by   simple    integration    and  h(x)    and  K(z)    are 
found   from   solving   LaPlace's   Equation: 

KW=A,Co2»K>x  +  A2sl*wAx 

(65) 

[•Mr  B,  Cos      A  i     +    &2^>u\A£ 

(66) 
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where  /*-  is  the  separation  parameter.   The  stress  function 
is  now  determined. 

~~£  yv  (67) 

Differentiating,  the  stresses  become 

<JX2  C,-X1(a,  CobhXx  +  Aj-SunkAxX^,  CosAi  +■  b2  Si*,  /i) 
CT£r    A^A.Co^Ax  +  AtSuvk>x)(&,Co*>£  +■  ^SL*  >i)     (68) 


The  normal  stresses  must  be  symmetric  in  x;  therefore, 

A2  =  0 
Without  loss  of  generality,  A^  =  1.   This  simplifies  the 
stress  expressions  somewhat  to  give, 

(Ti  =    A4Cosk>X  (5,  Cos  ^    +    £>z  Si.a.  >i) 
=  X>SUk>fc  (b,SUAt-  B^CosAi) 


(69) 


'** 
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We  require  that  the  shear  stress  be  zero  on  the  free  surface 
z  =  0.   To  accommodate  this,B_  is  set  equal  to  zero.   The 
stresses  simplify  further  to  the  form 

GV  -  M  6,  Cos»k  >>*  £**  ^  (70) 


'* 


*i 


>^&  iukM  6L**i 


Going  back  to  the  compatibility  equations,  others  can  now  be 
satisfied.   To  start  with, 

Q    -   G;  *  <Ty      +■  <Tft   ~    d  +  CTy  (71) 

Because   V  Q".  s  q   from  equation  (70),  equation  (56)  and 
(71)  give 


(72) 


Integrating, 


(T     r  X  FU)  +  6U)  (73) 
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Substituting  into  equation  (56) 


x  F"U)  +  G"^)  = 


o 


(74) 


The  only  way  that  this  can  be  satisfied  for  arbitrary  values 
of  x  and  z  is  for 


F"(*)  =  o 


and 


6"  (2.)  =o 


75) 


Then  upon  integration, 


Ct(*)  =  Cu£^C7    l76) 


and 


£7  -   %  tc**+cs)  f  (c^  4-  c.) 


(77) 


Since   \7  U,  -  O    from  equation  (70)  ;  then  equation  (56) 
gives 

}hG  .  JL  [c,  +  x(c4i^s)^  (A*.+  C7)  =o      (78) 


J?     ^L 


j 


This  is  satisfied  for  all  values  of  the  constants.   Since 
y  T\i  -  O         from  equation  (56) 


>X  ^i.         ^4 


c,  ■»-  * 


C4*  +  Cs)  +  (c«.44-Ci)/0 


(79) 
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This  equation  is  satisfied  if  C.    =  0;  therefore,  all  of  the 
field  equations  have  been  satisfied  and 


Cy    -     H*X  ^W^  C7 


(80) 


The  strains  in  the  plane  of  the  coating  are  now  of  the  form 


Ec 


Vt 


(81) 


C<T*  + 


Cu*  +  C7-T)cCt 


The  loading  boundary  conditions  are:  £^-  ^^^    cy  -  £y^  <tt  £--~Kc' 
This  gives 

(82) 


To  satisfy  these  requirements,  we  set 


At  JL       asJ.      cs-o 

2k0 


(83) 
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Cg  is  set  tc  zero  to  seek  a  solution  for  a  constant  ay. 
Since   Cnir-\)s  £y*>   ,  equations  (82)  simplify  to 

-  \)5  &c  £ys  =   C- 1  "  ^  ^i 

Solving  for  the  constants 


(84) 


Cn  =  £c  ^y*  Q-  ^0 


(85) 


The  last  constant  B,  is  found  by  setting  the  average  trans- 
verse stress  equal  to  zero  at  the  boundary  of  the  specimen. 
Assuming  the  specimen  to  be  W  units  wide,  0x:o  At  ^-   '  2-  • 


From  which, 


fe,=  "^c^ifgy 


(87) 
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The  stress  expressions  in  final  form  are: 


*  6-V) 


cHuTX 


se 


(88) 


a      and  t    are  neglected  as  being  approximately  zero 
compared  to  a  .   As  can  be  seen  from  the  expression,  if 
v   *  v  .  the  stress  state  becomes  uniaxial  in  the  limiting 
condition;  however,  if  there  is  a  Poisson's  Ratio  mismatch, 
and  there  normally  is,  a  transverse  stress  is  created  in  the 
coating  which  influences  to  some  degree   the  accuracy  of 
the  readings.   This  can  be  accounted  for  through  this  analy- 
sis by  formulating  an  expression  for  an  equivalent  Poisson's 
Ratio.   To  do  this  we  first  formulate  expressions  for  the 
average  strains  through  the  thickness  of  the  coating,  which 
will  be  what  the  fringes  will  be  measuring. 


(»-w) 
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An  equivalent  Poisson's  Ratio  is  defined  in  the  obvious  way. 

7  ** 

V    -      -      —  (90) 

Substituting   elastic   material  properties   referred    to   earlier 


,»«v       '  (91) 

UoV        +       -0Z.04         __ 

n-  >^ 


c^(ra 


A  plot  of  this  expression  was  made  assuming  that  W>>h  ;  for 

instance,  let  W  =  20h  .   A  plot  of  the  behavior  is  shown  in 

c  r 

Fig.    7. 

If    the    specimen  material   is    taken    into    the   plastic    range   where 

-^      -      0.500,    the   expression   for  Poisson's    Ratio  becomes 

n   _ V  4k<-' (92) 


For  the  same  parameters  of  W  =  20h  ,  the  plot  of  equivalent 
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Poisson's    Ratio   for    the   plastic   case    is   given    in   Fig.    8. 
In   both    cases    the    combined   specimen  has    at    the    free   boundary 
the   Poisson's    Ratio   of    the    coating;    however,    this    falls   off 
rapidly   as   one   moves    to    the   interior.      Within    three   or   four 
coating    thicknesses    of   the   edge    the   Poisson's    Ratio  has 
approached   a   constant   value,    but   it   is    different   than    the 
specimen  value    alone. 

EXPERIMENTAL    COMPARISON    OF    STRESS    CONCENTRATION 
IN    A   STRIP    WEAKENED   BY    A    CIRCULAR   HOLE 

To   test   the    validity   and   accuracy   of   utilizing    the 
photoelastic   method,    a   test  case   using    a    circular  hole    speci- 
men of    finite  width  was    selected.      R.    C.    J.    Rowland    (8) 
determined   the    stress    distribution    around    the    circular   hole 
and  with    the   known   theoretical   solution,    it  was   possible    to 
determine    the    accuracy   of    the   photoelastic  method. 
EXPERIMENTAL   PROCEDURE 

A  specimen  was   designed   to   accommodate    a   2.0    inch    diamete: 
hole   in   the    center.      To   determine    a   reasonable  width   and 
length  which  would   allow   a   uniaxial   field   to  exist  between 
the   center   hole    and   the   end   grips,    a    1/8    scale   PS-lC  model 
was    tested   under   various    tensile    load   conditions.      It  was 
determined   that   a  minimum  of    a   five- to-one    ratio   for    length- 
to-width   must  exist   to  eliminate    interference   of    fringe 
patterns.      Fig.    9    shows    the   geometry    selected. 
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CIRCULAR  HOLE  SPECIMEN 


40 
inches 


2    Inch 
Diamete: 


* 8<  inches" 


•Data    Points 


t   =    0.0916    in 


Figure   9 
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The    aluminum   specimen  was   placed   under   an    8000    pound 
load   in   the   MTS   machine,    and   compensator    readings  were    taken 
with   a  Model   0  32   Photolastic   Compensator    at   the   90°    position 
on  each   side   of    the   hole.      They  were    found   to  be    122    /    121 
counts,    or   an   average   of    121.5    counts.      The    strain-optic 
equation   from  Appendix  A   can   be   written   as 

£      -U    -      RrQ+OQM  (93) 

At  the  circular  boundary,   £.x  -  -  Vc  %.y 

4   -   F<r^  _  o.oVbZ^S    u/tn.  (94) 

^ '  "SI  " 

Making  the  correction  for  reinforcement, 

£*-  =  \.oH,*>ly    -  0.00**>4V     U/ui       (95) 

The    local  hoop   stress    is 

Howland's   results    for    this    case   gave    34,719    psi    at   the   edge 
of    the   hole.      This    indicates   a   difference   of    2.4%.      Later 
we  will   show   that   this    difference    is   reduced  by  more    than 
a   factor   of    three   by  extropolating    to    determine    the    surface 
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stress  from  readings  in  the  interior. 
OBLIQUE  INCIDENCE 

In  order  to  determine  stress  values  away  from  the  edge 
of  the  hole,  oblique  incidence  was  used.   Fig.  10  is  a 
schematic  showing  the  coordinate  system  for  rotation  about 
the  Y  axis.   The  governing  equations  become, 


K 


i\i 


gg 


(97) 


(98) 


—  i 


Figure  10.   Oblique  Incidence  Coordinates 
for  Rotation  about  the  y  axis. 
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The  strain  in  the  x1  direction  is  given  from  the  strain 
transformation  laws 


tl=  £x  Co^e  +  tx  Sc^e 


(99) 


For   plane    stress,    the    following   stress-strain    relationships 
hold: 


U  '-  -^i  foe  ^yc 


(100) 


Ql.  - 


Ec 


*C 


l-1> 


t,   L 


fc«  *  °c  t 


,v 


(101) 


CTvc - 


£c 


l--0cL 


ty  ■*■  l)o  £*. 


(102) 


Putting   these    together, 


*  l.-Oc 


(10  3) 


and  the   expression  of   strain   along    the   oblique   axis   becomes 


U--      6*  <U*"©  -    T>«      £fc,+  ty)SU*© 


-0, 


(104) 
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From  Fig.    10    it    can  be   seen    that      V\      —    *\  C-cs&  .      Putting 
all   cf    these  expressions    into  Equation    (98) ,    an   equation    for 
N  £Ni    in   terms   of   known  quantities   can   be  obtained. 

M«,=  T^f ! ]\lJ\--dt.Co!>te)+lL^~C0>&)l      (105) 


Solving  this  equation  simultaneously  with  equation  (97), 
t,    z    _Jl VNJw6-0c)Co»©  -  M^-Octe*1©) 


(106 


1  "(|+i)«.)Si^© 


INJ, 


L 


v  .1 

^(i-o«.)Cc^-m(ccS©-i>c;j    (107) 


For   rotation   about    the  other    axis, 


si 


frx 


fey      -fe.) 


(108) 


Fig.    11   shows    a   schematic   of    this   oblique    incidence 
arrangement.      The    strain    in   the    oblique    direction    turns    out 
to   be   of    the    same    form   as   before . 


t 


ty     =        ty     dc^Q     +      tl      SU-© 


(109) 
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Figure  11.   Oblique  Incidence  Coordinates  for 
Rotation  about  the  x  axis. 


Using  the  expresion  for  £.£  previously  derived  and  the 
expression  for  K  /  the  expression  for  the  oblique  fringe 
value  becomes 


N    -^ 


(l-Oc)  Cos6 


£y(Co^e-Oc)-^(i-^Cos1©) 


Solving  this  simultaneously  with  the  normal  incidence 
equation  (97) ,  the  strain  expressions  are: 


t~- 


0+iV)SbteL  J 
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1     6  +■&.)  s-ui©L 


- 

M(i-OcCoS1©)- MdK(i-i)c)Co<,eJ    (112 


Calibrating  the  Polariscope 

In  oblique  incidence  the  determination  of  the  angle 
can  be  tedious  because  of  changes  in  it  at  the  surface  due 
to  the  index  of  refraction  being  different  in  the  coating 
than  it  is  in  air.   The  light  ray  is  bent  as  shown  in 
Fig.  12.   It  is  more  reliable  to  obtain  the  value  experimen- 
tally.  For  instance,  equations  (106),  (107)  (111)  and  (112) 
can  be  written  as 


Ft  [aNsx,  -  CM 


£y  *    Ft 


CM 


_    A  N 


e* 


(113) 


(114) 


(115) 


(116) 


Subtracting  (114)  from  (113)  or  (115)  from  (116) , 

ty-U-~    MFt(c-6) 

Since  C  -  3  =  1  from  the  basic  equation  (A-l) , 
C  =  B  +  1 


(117) 
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Figure  12.   Definition  of  Angle  of  Refraction 
for  the  Polariscope 
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Solving    two   of    the    other   equations    simultaneously    for   A  and  B, 
j\    -         1/ i_— r-  (118) 


ttr£±)(ti&1-^9t) 


(119) 


Using    a   uniaxial    test    to   make    the    calibration    tests,    the 
expression   for   A,    B    and  C   become 

A-- -J- (120) 

N^  -   0  Nay 

[*>  = r  (121) 

(i  +  o)(Ney-  isle*) 

£   -      &  +■  I  (122) 


Several  readings  were  made  of  normal  incidence  and  both 

oblique  incidence  fringes  on  a  uniaxial  specimen  of  coating 

material  at  different  load  levels.   N„   and  I<L   were 

9x       8y 

plotted  versus  N  and  a  linear  regression  was  made  of  the  data 

Mex  -    o.Sid    ^O.KtN  (123) 

£_  -    0.9969 
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M   =  -2.166  +-  I-1&4  M  (124) 

R_  -  0.9997 


Substituting  into  equations  (120)  to  (122)  and  using 
"0  =  v<,=  0.3793,  the  constants  were  found  to  be 

A  =  1.372 
B  =  0.899 
C   =    1.899 

Generalizing   for   Any   Material 

In  order   to   use    the   equations    of   oblique   incidence    for 
measurements    on   materials    of    any  Poisson's    Ratio,    the   ex- 
pressions   for    the    constants   were    used    to    solve    for    the   overall 
incidence   angle    0. 

/\  - '        -      {'>lL-  (125) 

(  I  4-0c)     Swv>0 

Solving    this   equation,    0   =    31.85    .      The   other   two   expressions 

for  B    and   C  both   gave   9   =    31.76°.      The   good   agreement  gave 

confidence   in    the    numbers,    and   the    average   value   of    31.8° 

will  be  used  in  all  further  calculations. 

The  expressions  for  strain  in  terms  of  N  and  N„   become 
^  8y 
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For  rotations  about  the  x  axis, 


t*    ~~ 


l+OL 


(l.icoit.-h.Loi-L:o)t4  -  3.o^7(i-5)rJg.K 


V" 


.     Ft 


+  0 


(^.fcoii.-Z.froito)  ^  - 


^.ck©7 


(j-^M** 


PHOTCELASTIC    DATA   FOR   THE   HOLE 

Normal   and   oblique   compensator   readings   were    taken   at 
0.0625    inch   intervals    at    the   90      position   en   each   side   of 
the  hole    as    shown   in   Fig.    9    for   a    loading   of    12,000    lbs., 
or   a   far-field   stress   of    16,124   psi.      Six   different  compen- 
sator  readings   were  made  by    two  observers;    i.e.    each   observer 
made   three    readings,    and   these    are    listed   in   Table   VI.      The 
average    and   standard   deviation   are    also    listed. 
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TABLE    VI 
CIRCULAR  HOLE    COMPENSATOR   READINGS 
NORMAL   INCIDENCE 


X 

AVG. 

STD    DEV 

0.0625 

148 

151 

149 

149 

15  4 

14  8 

149.8 

2.32 

0.1250 

125 

130 

126 

124 

130 

125 

126.7 

2.67 

0.1875 

10  8 

110 

108 

108 

108 

10  7 

10  8.2 

.98 

0.2500 

94 

96 

97 

98 

97 

97 

96.5 

1.38 

OBLIQUE    INCIDENCE 


X 

AVG. 

STD    DEV 

0.0625 

126 

125 

124 

125 

12  8 

124 

125.3 

1.51 

0.1250 

102 

102 

100 

102 

10  4 

101 

101.8 

1.34 

0.1875 

85 

85 

84 

85 

87 

86 

85.3 

1.03 

0.2500 

74 

75 

75 

75 

76 

75 

75.0 

.63 

By  Chauvenet's    Criterion,    data  will  be   excluded  where 
probability   of   occurrence  would  be    less    than     *-    .      For   six 
readings    this    is    .0  833.      Assuming   a   Gaussian    distribution 
of   errors   on   the    readings,    this    represents    1.73a    away   from 
the   mean  before    the    datum  becomes    suspect.      Table   VII    shows 
the    upper   and   lower    limits    of    acceptable    data. 
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X 
0.0625 
0.1250 
0.1875 
0.2500 


TABLE    VII 
LIMITS    OF   ACCEPTABLE    DATA 


NORMAL    INCIDENCE 
UPPER   BOUND       LOWER   BOUND 


153.8 

131.3 

109^9 

98.9 


145.8 

122.0 

106.5 

94.1 


OBLIQUE    INCIDENCE 
UPPER    BOUND      LOWER   BOUND 


127.9 

104.1 
87.1 
76.1 


122.7 
99.5 
83.5 
73.9 


Data   falling  outside   of    these   ranges   were   eliminated   and 

new  averages   computed.      These   results    are    listed    in  Table  VIII 


TABLE    VIII 

SMOOTHED    CIRCULAR   HOLE    COMPENSATOR   DATA 

COMPENSATOR   READINGS 
Cn  C 


DATA    POINT 

(in.) 

0. 

.06  25 

0, 

.125 

0, 

.1875 

0, 

,25  0 

14  9.0 

126  . 7 

107.8 

97.0 


ex 

124.8 

101.8 

85.3 

75.0 


Extrapolating    to    the   Edge    of    the   Hole 

Because    it    is   difficult    to  measure    the    fringe    count 
right  at   the   edge   of    the    hole,    these    data   were   plotted 
and   extrapolated  back  to   that  point.      It  was    found    that 
log    log   values   of    the    compensator   readings   plotted    linearly 
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with  x.   The  best  fit  was  obtained  by  using  the  first  three 
points. 

Jbn  (i*g  C©*)  -  0.3390- azassx  (130) 


with   R  =    -.9996.      Using    the  value    for   x  =    0, 


cex  =  152'31 

For  this  Compensator,  47  counts  constituted  1  fringe. 

N8X=  3-24 


Similarly   for    the    normal   incidence    fringes, 

jL3  (^  Cn)  =  0.*sn-o.z.«** 


(131) 


with   R  =    .99996.      Using    the    value   for    x  =   0, 

CX7   =    176.81 
N 

and 

N      =3.76 

At   the   edge   of   the   hole    the    radial   stress,    a    ,  is    zero, 

which   requires    for    the    two   dimensional  plane    stress  geometry 
that 

tr    -    ~5    te  (132) 

Substituting  equations  (128)  and  (129)  into  equation  (132), 
a  relationship  is  found  between  N    and  N  for  all  ~d   at 
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the  hole  boundary. 

Nftx  -  a  84  99  nJ  <133> 


After  extrapolation  the  measured  values  of  N  =  3.76  and 

N.   =  3.2  4  check  out  to  within  1.3%  of  the  theoretical  values 

0  X 

from  equation  (133)  . 

Accounting  for  Poisson's  Ratio  Mismatch 

These  fringe  values  are  extrapolated  from  interior 
regions  where  "0  =  .3256;  however,  at  the  edge  of  the  hole, 

"0  =    .3793;  therefore,  they  must  be  adjusted  to  the  appro- 
priate Poisson's  Ratio  value.   This  is  done  by  knowing  that 
at  the  boundary  the  hoop  stress  is  given  by 


i  +  -0 


*    _^ <134) 


For   a  given   afl    existing    at   the   edge   of    the   hole,    the 
relation  between    the   photoelastic   readings    for   a   different 
Poisson's    Ratio   is    given    from    (134)     to   be 

-.  N*"  (135) 


1  +  O-L 


This,  of  course,  assumes  E   and  F   are  held  constant.   Thus 

s   _   e 

for  our  case,  \)x    =  .3793  and  VL=  .3256,  and 

N 
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The  actual  normal  incidence  fringe  value  at  the  edge  of 
the  hole  is  now 


M  ^  (l.O40b)(Vll)r  3.91 


To  find  the  corresponding  value  for  the  oblique 
incidence  fringe,  equation  (133)  is  used  in  conjunction 
with  equation  (135) . 


M 


Ox 


r  I.D405  Nex 


J-&=.yi<ta 


(136) 


-Or.^iSt 


From  this  the  value  of  N_   at  the  edge  of  the  hole  is 


MeK  =   (i.o4os)(a.i4)=  5.*7 


The  last  correction  to  the  data  is  for  reinforcement, 
which  for  these  materials  is  a  factor  of  1.0162  from 
equation  (39)  . 


N  =  3.9  8 


NQ   =  3.43 

ex 


Check,  against  Numerical  Solution 

A  check  can  now  be  made  against  the  numerical  solution, 
which  for  this  geometry  and  loading  is  given  by  a  finite 
element  solution  developed  by  Kaiser  (4)  to  be  a      =  51,678  psi 
Using  N  =  3.98,  the  extrapolated  photoelastic  measurements 
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give  from  equation  (134). 

This  value  is  within  0.3%  of  the  finite  element  value. 
This  is  surprisingly  good  agreement. 
Check  of  Stress  Distribution 

A  check  of  the  accuracy  of  the  photoelastic  work  was 
also  made  on  the  stress  distribution  away  from  the  hole  by 
means  of  a  finite  element  solution.   The  criteria  established 
to  determine  uniform  boundary  stress  distribution  was  uni- 
formity in  nodal  displacements  along  the  loaded  edge  as 
discussed  by  Segerlind  (9).   In  the  models  used  for  FINITE 
ELEMENT  ANALYSIS  (FEA)  in  this  report,  nodal  displacements 
were  uniform  to  within  0.1%,  and  the  resulting  stress  distri- 
bution was  uniform  axially  to  within  0.1%  at  the  panel  ends. 
In  all  cases  two-dimensional,  eight  noded,  isoparametric 
elements  were  used.   These  higher  order  elements  cannot  be 
loaded  in  an  "intuitive"  manner,  as  discussed  by  Zienkiewicz 
(10).   Fig.  13  shows  the  nodal  loading  required  to  obtain  a 
uniform  surface  load. 
Element  Meshes 

Two  meshes  were  developed  for  each  panel  analyzed. 
A  reasonable  effort  was  made  to  keep  element  corner  angles 
as  close  to  90°  as  possible  to  reduce  the  effect  of  element 
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FIGURE    13 
NODAL      LOADING       DIAGRAM 
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distortion  discussed  by  Hopkins  and  Gifford  (11) .   All 
meshes  modeled  a  quarter  of  the  actual  panel  by  using  the 
two  axes  of  symmetry  as  is  common  practice  in  FEA.   The 
step  from  course  to  fine  element  meshes  was  made  so  that 
each  element  in  the  course  mesh  was  subdivided  into  four 
smaller  elements  of  the  same  type.   Such  a  mesh  subdivision 
can  be  expected  to  give  monotonic  convergence  of  results, 
as  shown  by  Cook  (12) ,  and  allow  extrapolation  to  results 
of  an  infinitely  fine  mesh.   Figures  14  and  15  illustrate 
the  element  meshes  used  in  this  analysis. 
Using  ADINA 

Once  the  mesh  was  developed,  input  data  was  prepared 
in  accordance  with  the  ADINA  user's  manual  (13) .   This  same 
set  of  data  was  used  as  input  for  PSAPl,  the  plot  routine,  to 
check  for  errors  and  provide  a  graphical  display  of  the 
element  mesh.   After  preprocessing  by  PSAPl,  the  data  is 
entered  into  ADINA  for  analysis.   In  the  case  of  linear 
analysis,  two  types  of  stress  output  may  be  specified, 
nodal  point  or  Gauss  integration  point.   Nodal  point  output 
can  be  computed  for  up  to  eight  node  point  stresses  for 
each  element.   Since  2x2  Gauss  integration  was  used,  four 
Gauss  point  stresses  were  computed  for  each  element.   The 
2x2  Gauss  integration  is  recognized  as  the  most  efficient 
integration  order  for  this  type  of  analysis  [Ref.  10,  p.  2  84] 
The  linear  analysis  used  an  isotropic  linear  elastic  material 


65 


o 


2b  ELEMENTS   ( ISOPARAMETRIC) 

111  NODES 

IlJ2    DEGREES    OF    FREEDOM 

DIMENSIONS 


X 

W                   L 

RADIUS 

II                  1) 

0.2                5                25 

i 

1 

■i   i           !•                        H 
0.25  :    4.0625           20                        1 

i                !                       !                      ! 

X   = 


RADIUS 

w 


SYMMETRY 
AXES 


FIGURE  14 
COURSE    MESH    FOR   CIRCULAR   HOLES 
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A   SUBDIVISION    OF    COURSE    MESH 

112    ELEMENTS     (ISOPARAMETRIC) 
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FIGURE   15 
FINE    MESH    FOR    CIRCULAR    HOLES 
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model  (MODEL  "I"  in  Ref.  4)  which  required  input  of  tensile 
modulus  and  material  properties. 

For  static  analyses  ADINA  uses  a  time  function  method 
to  apply  loads  in  steps.   Linear  analysis  loading 
was  accomplished  in  a  single  step  to  a  nominal  value  of 
3,000  lbs  load.   The  stress  output  from  ADINA  is  a  listing 
of  nodal  or  Gauss  point  stresses  for  each  element.   Since 
the  only  area  of  interest  in  this  analysis  was  the  distri- 
bution of  stresses  along  the  reduced  cross-section,  no  large 
post-processing  program  was  developed  or  used.   All  final 
computations  using  ADINA  output  data  were  accomplished  on 
a  HEWLETT-PACKARD  9 8 30 A  calculator,  using  short  programs 
coded  in  BASIC.   If  more  extensive  stress  distribution 
information  were  desired,  some  form  of  automated  post- 
processing would  be  necessary  to  reduce  the  computational 
workload.   At  a  minimum,  nodal  stress  outputs  by  ADINA  must 
be  averaged  to  obtain  unique  values  of  stress  at  nodes  shared 
by  more  than  one  element. 
Richardson  Extrapolation 

The  use  of  course  and  fine  meshes  allows  extrapola- 
tion to  an  infinitely  fine  mesh  by  means  of  the  Richardson 
extrapolation  (14)  where: 


/A         /    \m 


G 


<£(kf)  -  <rF(kc) 


***f  "        c  -  K 


^ (137) 

c 
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where 

Gextrotp     =   extrapolated  solution 

(Tf  =   solution  obtained  with   h=h 

w  w  c 

<Te  =   solution  obtained  with   h=h_ 

F 

Wc  =  linear  dimension  of  course  element 
y\p  =    linear  dimension  of  fine  element 
^  =  2  (for  this  analysis) 
The  exponent  m  is  determined  by  the  order  of  the  discre- 
tization error  0(h  )  .   Since  h  represents  the  length  of  an 

2 
element,  the  element  area  is  represented  by  h  .   In  a  two 

dimensional  problem  such  as  this,  0(h  )  is  of  the  order  of 

2 
h  ,  the  area  of  an  element.   In  the  mesh  refinement  scheme 

used  Equation  (137)  can  be  rewritten 


Gwt^p  =    —  (138) 


thus 


(T  -         °P  _  4 


V^ro.^ 


4 


(139) 


Equation  (139)  then  becomes  the  relation  to  obtain  extrapola- 
ted stresses  from  coarse  and  fine  mesh  results  in  a  two  dimen- 
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sional  analysis.   Better  extrapolations  can  be  obtained  by 
using  three  or  more  defined  meshes,  but  the  computational 
effort  increases  significantly. 
Optimal  Stress  Locations  and  Local  Smoothing 

It  is  generally  accepted  that  the  most  accurate 
sampling  points  for  stresses  are  the  Gauss  integration 
points  within  the  element  [Ref.  10,  p.  281].   In  this  analy- 
sis, the  nodal  points  are  of  the  greatest  interest;  thus 
a  technique  of  local  smoothing  must  be  applied  to  the  inte- 
gration point  stresses  to  obtain  nodal  stresses  as  reported 
by  Hinton  and  Campbell  (15) .   The  nodal  values  obtained 
must  then  be  averaged  if  shared  by  two  or  more  elements. 
Circular  Holes  in  Linear  Material 

The  FEA  results  for  a  circular  hole  in  a  finite  width 
strip  were  used  to  validate  the  photoelastic  procedures  dis- 
cussed earlier.   The  results  of  Rowland  (8)  were  compared  to 
both  the  Gauss  point  smoothed  results  and  the  nodal  output 
results  in  Fig.  16.   The  stress  concentration  factor,  o/a^, 
is  referenced  to  the  far-field  stress.   The  smoothed  results 
give  the  best  match  to  the  results  of  Howland  at  the  edge 
of  the  hole,  and  the  only  significant  variation  between  the 
two  FEA  methods  occurs  within  the  first  0.25  inches  from  the 
edge.   In  order  to  obtain  the  0.25  inch  stress  value  for  the 
coarse  mesh,  in  the  Gauss  point  smoothed  result,  a  midside 
node  value  had  to  be  obtained  by  the  averaging  method. 
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TABLE    IX. 

A   =    0.25    FEA 

.    RESI 

JLTS    -    NODAL   OUTP1 

DISTANCE 

Geo     = 

4500 

psi 

FROM  HOLE,    in. 

aQ    ,psi 

ar    ,psi 

0.00 

14745.5 

329.8 

0.25 

9447.0 

2181.1 

0.50 

7232.6 

1599.3 

0.75 

6304.7 

1473.5 

1.00 

5771.7 

1139.6 

1.25 

5471.2 

891.6 

1.50 

5266.9 

671.0 

1.75 

5091.7 

489.1 

2.00 

4940.2 

354.5 

TABLE    X. 

A   =    0.25    FEA   RESI 

JLTS    -    GAUSS    OUTP' 

DISTANCE 

am   =    4500 

psi 

FROM   HOLE,    in. 

aQ    ,psi 

ar    ,psi 

0.00 

14422.8 

121.2 

0.25 

8721.2 

1927.1 

0.50 

7317.1 

16  39.7 

0.75 

6218.6 

1452.4 

1.00 

5779.6 

1151.5 

1.25 

5453.3 

891.0 

1.50 

5267.1 

677.3 

1.75 

5085.9 

486.9 

2.00 

4937.9 

343.8 
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The    linear   distribution   of    smoothed  stresses    along    the 
sides   of    the   element,     (15) ,    appears    to   produce   a   less    accurate 
result   in    this    area   of   extreme    stress   gradient,    when   compared 
to   ADINA's    nodal   output   result.      This    tendency  was   noted 
in   all   cases;    however,    the   peak    stress    values    from   smoothed 
results   consistently   gave   better   correlation   with   other    inves- 
tigators   ( 16)  . 

A   circular  hole  with   A=0.25    was    also   analyzed   and 
compared   to    the   photoelastic  experimental   results.      The 
stress   distribution    is   given   in   Table    IX   and   Table    X.      The 
data   for    the    first   two    locations    has   been   normalized   and 
entered   into   Table    XI    to   show   a    comparison  with   photoelastic 
results    similarly   normalized. 

TABLE    XI 
COMPARISON    OF   STRESS    DISTRIBUTION 

Finite  Element  Photoelastic 

x  ar/aoo  ax/a°°  ar/ooo  cg/a^, 


0.00 

-.0269 

3.2051 

-.116  7 

3.0780 

0.25 

.4282 

2.0993 

.4111 

2.1627 

Comparisons    of    the   photoelastic    data  with    the    finite 
element  show  agreement  within    3%    and   4%    on    the    average.      This 
is   not  as   good   as    the    correlation  made    at   the   edge   of    the 
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hole  where  we  required  a   =0.   If  a   is  allowed  to  be 

i       r  r 

whatever  the  readings  show;  then  the  correlation  is  in  the 
3%  or  4%  range;  however  as  can  be  seen  at  x  =  0 ,  the  finite 
element  is  not  exact  either,  especially  at  these  small  x 
values. 
Sensitivity  Studies 

At  the  edge  of  the  hole  where  "0    =  .379  3,  the  expres- 
sions for  strain  from  equations  (12  8)  and  (129)  become 

t     -I.  S198  X»cr*M  -  Z.4Z-9k  *  »°~*  ^**      (140) 


i     --   3.W*6*i°"*M  -2.4L%kio'  ^e* 


(141) 


Substituting  these  expressions  into  the  plane  stress  consti- 
tutive equations 


(142) 


<Jr  -   2. 08*51  xio4M    -  ZAS31  xio^Ne* 

<T   -   Z.«95^io4N    -  Z.4s£7x./o4nU*  (143) 


The  sensitivity  of  the  stress  is  given  by 


(144) 
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Applying  this  to  equations  (142)  and  (14  3) 

<d(Tr  -  Z-0891  *<oMnI    -    Z.4S67  *>i&4cl^     (145) 
JcT   =  Z.  £9  30  *  lO4^  JsJ    -   Z.A-SSIk  lo^d^e^      (146) 


Selecting  the  maximum  deviation  of  N  and  N   as  1.73a,  where 

d  X 

a    is    the  standard  deviation   of    the   smoothed   data 


dM  -  (i.i5)(t  1.22)  =  ±  z-n 


Then  da  can  have  two  values  depending  upon  how  the  plus  and 
minus  signs  couple  up.  The  worst  of  these  are: 

<4<Tr  -     -    i 6,9 7    p^ 
<d(Te  =    -    Zo""l    f^ 

These   figures    represent  the  maximum  errors    that  would  be 
expected   due    to  misreading    the   fringe  pattern,    and   they   are 
approximately   4%    on   a      and  explains   how  non-zero   results 

for   a      arise. 

r 
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CONCLUSIONS  AND  RESULTS 
From  very  careful  experimental  procedures  where  no 
handbook  values  for  material  properties  were  used  but  all 
were  carefully  determined,  photoelastic  measurements  checked 
theoretical  results  within  4%.   This  error  included  all  errors 
in  material  properties,  in  operator  readings  and  in  equip- 
ment calibration.   Without  very  careful  and  detailed  measure- 
ments on  every  parameter,  photoelastic  results  could  easily 
be  in  error  by  15%  to  20%. 

With  this  level  of  accuracy  established  by  controlled 
laboratory  procedures  and  calibrations  of  equipment  in  place, 
we  will  now  proceed  to  study  the  elastic  and  plastic  behavior 
of  notches  where  fatigue  failure  originate. 


76 


REFERENCES 


1.  Frank   E.    Stenstrom,    "Photoelastic  Study    of   Elastic    and 
Plastic   Stress    Fields    in    the   Vicinity   of    a  Notch", 
Master's    Thesis,    Naval  Postgraduate   School,    Monterey, 
California,    December    19  80. 

2.  Amineo-Tuckerman   Optical  Strain   Gage   System  Instructions 
No.    750,    p . 16 ,    American   Instrument  Company,    Inc.    1958. 

3.  ASTM    19  80    Annual  Book    of   Standards,    Part    10,    Standard 
E8-77a. 

4.  Michael   J.    Kaiser,    "An  Elastic-Plastic   Finite   Element 
Analysis    of   Notched  Aluminum  Panels",    Master's    Thesis, 
Naval  Postgraduate    School,    Monterey,    California, 
March    19  81. 

5.  Micro-Measurement,    Transverse   Sensitivity  Errors    TN-137, 
undated,    Vishay   Intertechnology ,    Inc. 

6.  Walter   Ramberg    and  William   R.    Osgood,    "Description   of 
Stress-Strain   Curves   by   Three   Parameters",      NACA  TN   902, 
July    1943. 

7.  J.    W.    Dally   and  I.    Alfirevich,    "Application   of   3ire- 
fringent  Coatings    to   Glass-fiber-reinforced  Plastics" , 
Experimental  Mechanics,    Vol  9,    p.    97-102,    1969. 

8.  Howland,  R.C.J. ,  "Stress  in  a  Strip  of  Finite  Width, 
Which  is  Weakened  by  a  Circular  Hole,"  Philosophical 
Transactions  of  the  Royal  Society,  Series  A,  V.  229, 
pp.    49-86,    London,    19  30. 

9.  Segerlind,    L.J.,    Applied  Finite  Element  Analysis,    p.    251, 
Wiley,    1976. 

10.  Zienkiewicz,    O.C.,    The    Finite  Element   Method,    3d  Ed., 
McGraw-Hill   Ltd,    1977. 

11.  David  W.    Taylor   Naval   Ship   Research    and   Development 
Center   Report,    Accuracy   Loss    in    Distorted  Isoparametric 
Finite   Elements,    by   d/a.    Hopkins    and   L.    N.    Gifford, 
September    19  7  8 . 

12.  Cook,    R.D.,    Concepts    and  Applications    of   Finite   Element 
Analysis,    pp.    87-92,    Wiley,    1974. 


77 


13.  Massachusetts    Institute   of   Technology   Report    824  48-1, 
A   Finite   Element  Program    for   Automatic    Dynamic   Incre- 
mental  Nonlinear   Analysis    (ADINA)~    by  K.    Bathe, 
September    19  75    (revised  May    19  76) . 

14.  Crandall,    S.H.,    Engineering   Analysis    -   A  Survey    of 
Numerical  Procedures,    pp.    171-173,    McGraw-Hill,    l9~5  6 . 

15.  Hinton,    E.,    and  Campbell,    J.S . ,    "Local    and   Global 
Smoothing   of    Discontinuous   Finite  Element   Functions 
Using    a   Least  Squares   Method,"    International  Journal 
for  Numerical   Methods    in  Engineering,    V.    8,    pp.    461-480 
1974. 

16.  Peterson,    R.E.,    Stress   Concentration   Factors,    pp.    20-36 
and   111-150,    Wiley,    1974. 


78 


APPENDIX  A 


Equations  of  Photoelas ticity  Coatings 


+  t^    kL 


(A-i) 


cr,  -^i.  c  1fr  ~  (A"2) 


(A-4) 


^  "   ll  +  Oci  (A-5) 
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